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Abstract 

The rare decay B + — > DfK*° can occur only via annihilation diagrams in the 
Standard Model. We calculate the branching ratio in perturbative QCD approach 
based on kr factorization theorem. We found that the branching ratio of this decay is 
about of order 10~ 8 , which may be sensitive to new physics. 



1 Introduction 

Rare B decays are very important in particle physics, because they are important windows in 
testing the standard model (SM) and they are sensitive to new physics. As a rather simple 
method, factorization approach is accepted, because it explained many decay branching 
ratios successfully pQ. Recently, some efforts have been made to improve their theoretical 
application. One of these methods is the perturbative QCD approach (PQCD), in which 

many rare branching ratios such as B — > [2], B — > tttt [3j, B — ► up j3] were predicted. 
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In recent calculations, B — > P4*" 1 i^W" 1 ", P + — > K have been analyzed in the 

PQCD approach Ej, leaving P + — >■ Z?g if not calculated. In decay B + — > DfK , 
none of quarks in the final states is the same as one of the B meson. This decay is a pure 
annihilation type decay, which is described as B meson annihilating into vacuum and Df 
and K produced from vacuum afterwards. Some information about PQCD picture can be 
get in [HJ. 

In PQCD, the decay amplitude is usually separated into soft($), hard(if) and harder(C) 
dynamics by different scales, the factorization theorem allows us to write the decay amplitude 
as convolution, 

Amplitude- Jd 4 hd 4 k 2 d 4 k 3 r Ti[C{t)^ B {k l )^ Ds {k 2 )^K*{h)H{k u k 2 ,k^t)\. (1) 

In the function, ki(i = 1,2,3) are momenta of light quarks in each meson. C(t) is Wilson 
coefficient which comes from the QCD radiative corrections to the four quark operators. $m 
is the wave function which describes the inner information of meson M. H describes the 
four-quark operator and the quark pair from the sea connected by a hard gluon whose scale 
is at the order of Mb, so the H can be perturbatively calculated. The hard part H is channel 
dependent, while $a/ is independent of the specific processes. 

Some analytic formulas for the decay amplitudes of B + — > D+K decays will be given 
in the next section. In section |31 we give the numerical results and discussion. Finally, we 
conclude this study in section HI 



2 B + -> DfK*° amplitudes 

For simplicity, we consider B meson at rest. In the light-cone coordinate, the B meson 
momentum Pi, momentum P 2 and K momentum P3 are taken to be: 

P 1 = -|(i,l,0 T ), P 2 = -|(l-r 2 ,r 2 ,0 T ), P 3 = -|(r 2 ,l-r 2 ,0 T ), (2) 



where r 2 = M d +/Mb and r% = Mj^o/Mg. Putting the light (anti-)quark momenta in B~ 
and K mesons as ki, k 2 , and k$, respectively, we can choose 



kx = (x!P+, 0, k 1T ), k 2 = (x 2 P+, 0, k 2T ), k 3 = (0, X3P3-, k 3T ). (3) 



— *0 

The K meson's longitudinal polarization vector e and transverse polarization vector 6t are 



given by : 



Mfj r, On Mb , s . . 

( ^v^y rli - rl ° T) - e ^vs^ (o '°' 1t) - (4) 



Then, integration over k 1 , k 2 , and fc^" in eq.Q leads to: 



Amplitude ~ / dxidx 2 dx 3 b 1 dbib 2 db2b 3 db 3 

Tr [C(t)<S> B ( Xl , &i)$D a (x 2 , 6 2 , e)$jc.(x 3> & 3 )#(si, 64, e, t)^) e - s(t) ] , (5) 

where 6j is the conjugate space coordinate of kiT, and i is the largest energy scale in H, as the 
function in terms of Xi and b{. The large logarithms (lnm^/t) coming from QCD radiative 
corrections to four quark operators are included in the Wilson coefficients C(t). The large 
double logarithms (In 2 Xj) on the longitudinal direction are summed by the threshold resum- 
mation [9 , and they lead to a jet function S t (xi) which smears the end-point singularities 
on Xi. The last term, e contains two kinds of logarithms. One of the large logarithms is 
due to the renormalization of ultra-violet divergence In tb, the other is double logarithm In 2 b 
from the overlap of collinear and soft gluon corrections. This Sudakov form factor suppresses 
the soft dynamics effectively [10j . Thus it makes perturbative calculation of the hard part 
H reliable. 

The heavy B and D s meson wave functions are restricted by heavy quark symmetry. In 
the heavy quark limit, we may use only one independent distribution amplitude for each of 
them 

$ M (x, b) = [( A75) + M 75 ] 4> M {x, b), (6) 

where M — B, D s . For the light K* meson, only the longitudinal wave function for outgoing 
state is relevant, which is written as: 

$K*{x3,b 3 ) = -jJ\M K * fa<t>K*{x3,b 3 )+ 4l 54(^,M + M^%(x 3 ,i3)]. (7) 

Unlike the heavy mesons, there are three different distribution amplitudes. 

In the decay B + — > D^K , the effective Hamiltonian at the scale lower than My/ is the 
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same as B + — > DfK decay 



H, 



cff 



Gp 

71 



K* 6 K d [Gi(//)0 1 (//) + C 2 (//)0 2 (/i)] 



Oi = {h^P L d){c^P L u), 2 = (b^P L u)(c^P L d), 



(9) 



where the projection operator is defined by Pl = (1 — 7s)/2. V* b V C( i are the products of the 
CKM matrix elements, and Cx^il 1 ) are the Wilson coefficients. According to the effective 
Hamiltonian, the lowest order diagrams of B + — > D+K are drawn in FigHJ As stated 
above, this decay only has annihilation diagrams. 





d s d s 




(<0 (d) 



Figure 1: Diagrams for B + — > D+K decay. The factorizable diagrams (a),(b) contribute 
to F a , and nonfactorizable (c), (d) do to M a . 

After perturbative QCD calculations, we get the factorizable contribution F a , and the 
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nonfactorizable contribution M a illustrated by Fig.l(a)(b) and Fig.l(c)(d), respectively. 



pi poo 

F a = 1§-kC f M 2 b / dx 2 dx 3 / b 2 db 2 b 3 db 3 4>D s {x 2 ,b 2 ) 
Jo Jo 



o jo 

{(1 -x 3 - 3r 2 + 2x 3 rl)(f) K *(x 3 ,b 3 ) +r 2 (1 - 2x 3 )r 3 <f) t K *(x 3 ,b 3 ) 
- r a (3 - 2x 3 )r 3 <fi s K « (x 3 , b 3 )}E f (t 1 a )h a (x 2 , x 3 , b 2 , b 3 ) 

~ {(! - rl)x 2 (j) K * (2:3,63) 

- 2r 2 (l + x 2 )r 3 (j) s K , {x 3 , b 3 )}E f {t 2 a )h a {\ - x 3 , 1 - x 2 , 6 3 , &a)l , (10) 



1 r 1 r 00 

M a = —==QAnC F M 2 B I dxidx 2 dx 3 / M61 6 2 d& 2 (f)B{xi,b 1 )(j) Ds (x 2 ,b 2 ) 
\/2N c Jo Jo 

x [x 2 (l - 2rl)<p K *(x 3 M) +r 2 (x 2 + x 3 - 1) r 3 0^*(x 3 , b 2 ) 

- r 2 (x 2 -x 3 + l)r 3 (p s Kf (x 3} b 2 )}E' m (t] n )h^\x u x 2 ,x 3 MM) 
- {(-(^2 + 2x 3 - 2)r\ + x 3 ) (p K *(x 3 ,b 2 ) +r 2 (l-x 2 - x 3 ) r 3 (p l K , (x 3 , b 2 ) 
-r 2 (3 + x 2 - x 3 )r 3 (p s K ,(x 3 ,b 2 )}E' m (t 2 m )h {2) (x 1 ,x 2 ,x 3 ,b 1 ,b 



(11) 



In our work, r 2 , r\ and x\ in numerators are neglected. In the functions, C F = 4/3 is the 
group factor of SU(3) C gauge group, and the functions Ef, t^' 2 , h a are given in the appendix 
of 0. 

Comparing with the previous calculated B° — > Di*^ K + decay [3], we find that the 
leading twist contribution, which is proportional to 0^-, is almost the same. However the sub- 
leading twist contribution, which is proportional to r 2 r3 in (fTU|) and (fTTj) . is quite different. 

The total decay width for B + — > DfK decay is given as 

T(B+ -> DfK*°) = -§-^(1 - r 2 2 )\V: b V cd {fBF a + M a )\ 2 . (12) 

The decay width for CP conjugated model, B —> D s K , is the same value as B + — > 
DfK , just replace V* b V c d with V u bV* d . Since there is only one kind of CKM phase involved 
in this decay, there is no CP violation in the standard model. 
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3 Numerical Results 



We use the same B and D s meson wave functions as before j3] 



0b(x, b) = Nbx 2 (1 — x) 2 exp 



M 2 x 2 



(13) 



<t>D B {x,b) 



■J Da x{l - x){l + a Ds (l - 2x)}exp 



^d s &) 2 



The K* meson's distribution amplitudes are given by light cone QCD sum rules 



:6x(l - x) {l + 0.57 ■ (1 - 2a?) + 0.07 • C 2 3/2 (l - 2x)| 
: |o.3(l - 2x)[3(l - 2x) 2 + 10(1 - 2x) - 1] + 1.68C 4 1/2 (1 - 2s) 



(14) 



(15) 



J K 



*(x 



2^2N C 

fT 

JK* 

2^2N C 

+0.06(1 - 2x) 2 [5(l - 2x) 2 - 3] + 0.36[1 - 2(1 - 2x)(\ + ln(l - x)}} , (16) 



A'* 



{3(1 - 2a;) [1 + 0.2(1 - 2x) + 0.6(10s 2 - 10a; + 1)] 



2^2N C 

-0.12x(l -x) + 0.36[1 - 6x - 21n(l - ar)]} 



(17) 



with the Gegenbauer polynomials, 



c 2 1/2 (0 = ^(3e 2 -i), ci /2 (o = \m 4 -m 2 + s). (is) 

In addition, we use the following input parameters, for meson decay constants and the CKM 
matrix elements and the lifetime of B + |13j . 

f B = 190 MeV, f Ds = 220 MeV, /J£? = 200 MeV, (19) 
I Kb I = 3.6 x 10" 3 , |Kd| = 0.224, r B+ = 1.67 x 10" 12 s. (20) 

The branching ratio obtained from the analytic formulas may be sensitive to many pa- 
rameters especially those in the meson wave functions |5J Ej. Similar to the B — > D^K 
decays |3J El > we found that the branching ratio is not sensitive to the parameter change of 
the K* meson wave function, but they do be sensitive to the heavy B and D s meson wave 
function parameters. For illustration of the uncertainties of the branching ratios, we choose 
the following B and D s meson wave function parameters 

0.35GeV < uj b < 0.45GeV, (21) 
0.21GeV < uj Ds < 0.30GeV, (22) 
0.21GeV < a Ds < 0.30GeV. (23) 
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Table 1: Decay branching ratios calculated in PQCD approach using the same parameters 
for B and D s wave functions, and fo 8 = J'd* — 220MeV. 
Decay channel D+K° D* + K° DfK*° 
Br ( 10~ 8 ) 1.5 ±0.2 3.7 ±0.5 1.8 ±0.3 

Using the range stated above, the branching ratio normalized by the decay constants and 
the CKM matrix elements results in : 

Br(B+ - DtK*°) = (1.8 ± 0.3) x lO" 8 ( ^ ' ( Vcdl ) ' . (24) 

v ; v ; V!90 MeV-240 MeVy VO-0036 • 0.224; v ' 

Considering the uncertainty of fr> s and the CKM matrix elements, the branching ratio of 
the B + — > DfK ° decay is at the order of 10~ 8 . This is still far from the current experimental 
upper limit [T3] . 

Br(fi + -> D+K*°) < 5 x 10~ 4 . (25) 

For comparison, we list the branching ratios of B + — > Ds* K and B + — > DfK 
calculated in PQCD approach, using the same parameters for B and wave functions. 
For simplicity, we set fo s = /d* = 220MeV. From the above result, we can see that the 
branching ratio of B + — ► D+K is a little larger than that of B + — > DfK From 
equation (jlOlllj) and formulas in ref. 0, we find the coefficients of sub-leading twist become 
negative, but they are all proportional to r 2 r 3 which is a little bit small. Because of the 
suppression of CKM matrix elements, the branching ratio of B + — ► DjK is much smaller 
than that of the neutral decay B° — > D~K . 

Since pure annihilation type B decays are suppressed comparing to spectator diagram 
decays, the soft final state interactions may be important In our case, B + meson can 
decay into D*° and tt + , the secondary particles then exchanging a K + , scatter into Z?+, 
K through final state interaction afterwards. This picture is depicted in Fig|21 which is 
difficult to calculate accurately, since final state interaction is purely non-perturbative [T4] . 
In ref.[5j, the results from PQCD approach for B° — > D~K + decay were consistent with the 
experiment, which shows that the soft final state interaction may not be important. 
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K*°(sd) 




Dt(cs) 

Figure 2: Typical diagram for picture of final state interaction. 

4 Conclusion 

In hadronic two-body B meson decays, the energy release is larger than 1 GeV. The final 
state mesons are moving very fast. The soft final states interaction may not be important 
in the two-body B meson decays. The PQCD approach based on factorization theorem 
is applicable to the calculation of B meson decays. 

In this work, we calculate the B + — > D^K decay in the PQCD approach. Since neither 
of quarks in B + appeared in the final state mesons, this process occurs only via annihilation 
type diagrams. From our PQCD work, the branching ratio of B + — > DjK is small in SM, 
which is around 10~ 8 . This branching ratio will be measured in the LHC-B in future. This 
small branching ratio predicted in the SM, makes the channel sensitive to any new physics 
contribution. 
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